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If the damping of a simple harmonic oscillator from a ther- 
mally random force is sufficiently strong, then the oscillator 
may become unstable. For a photon oscillator (radiatively 
damped by electric dipole moments), the instability leads to 
a low temperature Hepp-Lieb-Preparata super-radiant phase 
transition. The stable oscillator regime is described by the 
free energy of the conventional Casimir effect. The unstable 
(strongly damped) oscillator has a free energy corresponding 
to Dicke super-radiance. 

PACS: 78.60.Kn, 78.60.Fi, 78.70.-g 



I. INTRODUCTION 

In many theoretical treatments of condensed matter 
systems, the Hamiltonian consists of a sum of kinetic en- 
ergy and Coulomb potential energy. In reality, the trans- 
verse electromagnetic (radiation) fields are also impor- 
tant. If the transverse electromagnetic field is included 
in the description of condensed matter, then phenomena 
such as the Casimir effect jy], the Lamb shift |2j-[|] and 
the long ranged inter-molecular interaction forces 
can be explained employing standard low order quantum 
electrodynamic perturbation theory. If the coupling of 
the transverse electromagnetic field with matter is strong 
enough to produce an instability in the photon oscilla- 
tors, then a transition occurs from a normal radiative 
phase to a super-radiant phase. The Dicke model || has 
been studied as a prime example of this phase transition. 
In the Dicke model, many two level molecules interact 
with a single photon mode of the radiation field. The 
Dicke model has been extensively studied by many work- 
ers [p~o| p~3|l , each adopting a somewhat different method 
to explain the physical properties of the system. Hepp 
and Lieb [fl()| , analyzed the thermodynamic properties of 
the Dicke maser model and elucidated the mathematical 
structure of the super-radiant phase transition. 

In subsequent works of Wang and Hoe |T1J and Hepp 
and Lieb p3], the coherent states of Glauber [O] were 
employed for multi-photon oscillators in the Dicke model 
and the resulting phase transition has been investigated 
in great detail in these papers. The existence of such 
phase transitions in condensed matter systems has been 
investigated in a series of papers by Preparata and 
coworkers |15-2(]]. Preparata has proposed several ex- 



amples of super-radiant phases in a recent review p8| . 
Water || |of is a particularly interesting example. 

The purpose of this work is to present a very di- 
rect method for computing the thermodynamic region 
wherein super-radiant phase exists. When a condensed 
matter system interacts with the electromagnetic field, 
the field degrees of freedom may be described by a renor- 
malized thermodynamic potential. For small oscillator 
displacements (in the normal radiation phase) the effec- 
tive potential is parabolic. As the thermodynamic pa- 
rameters change, the effective oscillator frequency may 
turn imaginary. The resulting oscillator instability is 
the signature of the phase transition into a super-radiant 
state. 

In Sec. II, the statistical thermodynamic method for 
calculating the effective potential of a photon oscillator 
mode will be discussed. The conditions for the stability 
of the oscillator will be exhibited. In Sec. Ill, the dy- 
namics in real time of the oscillator will be discussed in 
terms of the photon mode propagator. The conditions 
for dynamic real time stability are shown to be the same 
as for thermodynamic stability. In Sec. IV, the example 
of a mesoscopic object placed in an optical cavity is con- 
sidered. If the mesoscopic object has a sufficiently high 
polarizability, then a cavity mode can become unstable. 
In the Sec. V, the nature of the non-linear forces which 
stabilize the oscillator mode in the super-radiant phase 
are explored. The shift in the equilibrium position re- 
quired to stabilize the oscillator in the low temperature 
super-radiant phase is exhibited. In the concluding Sec. 
VI, the notion of super-radiance with electric dipole cou- 
pling is related to the notion of coherent ferro-electricity. 
A material grain with a high electric polarizability, but 
otherwise without a net dipole moment, can be induced 
into a ferro-electric phase by an interaction with trans- 
verse electric fields. The transverse fields may induce 
a coherence in the electronic states then yielding a net 
electric dipole moment. 



II. STATISTICAL THERMODYNAMICS 

Let us consider a simple harmonic oscillator with mo- 
mentum P, displacement Q and frequency ujoc. The os- 
cillator is coupled to a thermal bath described by the 
Hamiltonian H . The free energy of the thermal bath 
J-q is then determined by 
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e -F /k B T = Trbath e -H /kB T_ 



(1) 



The coupling between the thermal bath and the oscillator 
is described by a force /. The total Hamiltonian is then 



H = \(P 2 +u J l Q 2 )+H-fQ. 



(2) 



which determines the total free energy T via 

e -r/k B T = Tr e -H/k B T (3) 

If Q and P are classical to a sufficient degree of accu- 
racy, then the the total trace may be approximated by 



Tr(...) 



(dPdQ\ 



(4) 



From Eqs.(l-4) one obtains 

p -(r-^ )/k B T 



( dPd Q \ e -P 2 /2k B T e -G(Q,T)/k B T r§\ 

\ 2-kU J 



i.e. 



2-Kh 2 



dQe 



-g(Q,T)/k B T 



(6) 



where 



G(Q,T) = 



k B Thi 



Tr bath e -( g -/Q)/ fc ^ 
Tr bath e -ff/W 



(7) 



Finally, if the force / is also classical to a sufficient degree 
of accuracy, then 



f Tr bath e -iH-fQ)/k B T 

-H/k B T 



JQ/k B T 



(8) 



I Tr bath e 
yielding 

Gclassical(Q, T) = (^Q 2 /2) - k B T In (e^/ kBT ) . (9) 

The "effective" thermally induced potential of the os- 
cillator is given by 

G(Q, T) = a4(Q 2 /2) - k B T In ( e fQ/^T\ j (10) 



where (without loss of generality) one may assume (for 
Q = 0) that (f) = 0. For small values of the oscillator 
displacement, Eq.(10) implies 



WS(Q,T)/Q 2 ) = (u; 2 /2), 



(11) 



- 2 = -L-(/ 2 )AbT. 



(12) 



The frequency shift (from lo^ downward to u>o) of the 
oscillator due to the interaction of the oscillator with the 
thermal bath is given by 



(13) 



The oscillator is "stable" if the shifted frequency is real; 
i.e. ujq > or equivalently 



fcfl^L > (J 2 ) (stable). 



(14) 



The oscillator is "unstable" if the shifted frequency is 
imaginary; i.e. ujq < or equivalently 



kBTuo' 2 ^ <{f) (unstable). 



(15) 



One may appreciate the dynamical significance of the 
shifted frequency by considering the spectral density 
Sff(uj) of the force fluctuations ]2^-|25|] defined by 



(/(*)/(*')> = 
For equal times 

(/ 2 ) 



S f f(u>) cosuj(t — t)du>. 



2 / S f f{uj)duj 
Jo 



(16) 



(17) 



The classical fluctuation response theorem asserts that 
'k B T 



5fte T(w + £(T 



(18) 



where r(£) for £ > is the dynamic damping coef- 
ficient for the oscillator. r(£) will be discussed in detail 
in Sec. Ill below. Here wc note that Eqs.(13), (17) and 
(18) imply that 



W_ - LO n = | - 
7T 



3?e T(uj + i0 + )dco. 



(19) 



If the integrated (over frequency) dissipative damping 
on the right hand side of Eq.(19) is sufficiently large, the 
the oscillator becomes unstable. The oscillator will be 
stable for 



<*C > I - 

7T 



die T(lu + iO+)duj (stable). (20) 



Eq.(20) has been here deduced by statistical thermody- 
namic reasoning. However, it is perhaps more simple to 
deduce Eq.(20) directly from the simple dynamical equa- 
tions of motion for the oscillator. 



where, for a classical force /, the shifted frequency c^o of 
the oscillator is given by 



2 



III. REAL TIME OSCILLATOR DYNAMICS 

Let us now consider a simple harmonic oscillator with a 
canonical momentum P, displacement Q and frequency 
Woo. The oscillator is driven by a small time varying 
external force Sf ex t{t). The total Hamiltonian, including 
the external force is then 



Ti-totai = H — Sf ext (t)Q 



(21) 



where H is defined in Eq.(2). In a phenomenological 
linear response theory, the small mean deviation from 
thermal equilibrium in the oscillator coordinate obeys the 
equation of motion 



SQ(t) + Lo^SQit) = Sf(t) + Sf ext (t) 



(22) 



where 5f(t) represents the induced change in the fric- 
tional damping force exerted on the coordinate by the 
thermal bath. 

For simple engineering purposes, one often considers 
that the frictional damping force has a simple local (in 
time) form 8f(t) ~ ~-f8Q(t). More generally, the damp- 
ing force will be non-local in time; i.e. 



*/(*) 



g(s)6Q(t-s)ds. 



(23) 



For complex frequency with 3m £ > 0, if 

Sfextit) =ne{Sf ext (0)e- iit }, (24) 

then 

SQ(t) = 5Re {D(()5f ext (0)e-^} (25) 

defines the dynamical oscillator propagator -D(£). From 
Eqs.(23) and (25) 



«/(*) = Re {H(QD{Q8f ext (0)e-« t } 



n(C)=iC/ Q{a)e*'da. = »cr(c). 



where 



From Eqs.(22), (25) and (26), one finds that 
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D(Q = 



w ^-c 2 -n(C) 



(26) 



(27) 



(28) 



i.e. n(C) is the "self energy" of the harmonic oscillator 
propagator. 

From the retarded (causal) nature |26| of Eqs.(23) and 
(27) follows the dispersion relation 



n(C) 



lo 3to Tl(u) + i0 + ) 

(w 2 - C 2 ) 



duo (29) 



Equivalently, in terms of the frequency dependent damp- 
ing coefficient T(lo + i0 + ), 



3?m n{u + iO+) = lo 5Re T(lo + i0+) (30) 
Eq.(29) reads 

As C — > 0, Eq.(25) reads as the thermal (static) re- 
sponse function 

dQ 



D(0) 



1 

dfext J x w 
2 _ , .2 



where 

wg = <4,-n(o). 

Equivalently, employing Eqs.(31) and (33) 



Wn = I - 

7T 



5Re r(w + i0 + )du 



(32) 
(33) 
(34) 



as previously derived in Eq.(19) from purely thermody- 
namic reasoning. 

If the dynamical oscillator propagator is presumed to 
be retarded, it too must satisfy a dispersion relation 



D(C) = 



lo Sto D(lo + i0 + ) 



cLo. 



(35) 



If the oscillator is stable, then causality is obeyed with the 
dissipative condition Sm D(lo + i0 + ) > 0. The response 
function -D(C) is then analytic in the upper half 3to ( > 
of the complex frequency plane. 

On the other hand if the oscillator is unstable, then at 
least one pole at £ = iw appears in the upper-half plane 
and the dispersion relation Eq.(35) for the dynamical os- 
cillator propagator is no longer valid. The occurrence of 
an instability in the system can be viewed as inducing 
the non-analyticity of the dynamical oscillator propaga- 
tor given in Eq. (28) . To locate an unstable oscillator pole 
in the upper half plane at C — we look for a zero in 
the denominator of Eq.(28); i.e. 

lo^ - {ivof - IL(iw) = 0. (36) 

From Eqs.(31) and (36) we seek a solution for w > of 
the equation 



w = | — 

7T 



w 2 Ke T(oj + i0+) 



du>. 



(37) 



/ (uj 2 + w 2 ) 

As a function of imaginary frequency in the range < 
w < oo, it is evident that the right hand side of Eq.(36) 
is smoothly and monotonically decreasing from 11(0) to 
zero, while the left hand side of Eq.(36) is smoothly and 
monotonically increasing from lo^ to infinity. If 11(0) > 
LO^, then the two smooth curves intersect at precisely one 
value of w. It is easily proved that £ = izu is an isolated 
zero in the denominator of Eq.(28). Thus D(() has an 
isolated pole in the upper-half plane. The condition for 
this pole to exist, i.e. 



CtC < I - 

7T 



3?e T(lo + i,0 + )du), 



(38) 



is the same as the existence of the oscillator instability. 
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IV. OPTICAL CAVITY INSTABILITY 

Consider an electromagnetic mode in an empty optical 
cavity p7|-|32[ . The electric field of the mode Eoo (r) may 
be related to the oscillator displacement Q of the mode 
by equating the mean energy stored in the electric field 
to the potential energy of the oscillator; i.e. 



1) [ | Eoo (r)| 2 d 3 r=(i)^Q 2 . 

57r / J cavity \ z / 



(39) 



In order to implement Eq.(39), one may choose 

Eoo(r) = ljooQ eoo(r) (40) 
where e oc (r) is normalized according to 



f | eoo (r)| 2 d 3 r=l. 

J cavity 



(41) 



Let us further suppose that cavity has been designed so 
that the mode electric field is localized in a spatial region 
of length scale L about the origin. Then 



(42) 



where n is a unit vector; (n • n) = 1 . 

Let us finally suppose that a mesoscopic object is 
placed within the mode neighborhood. The mesoscopic 
object interacts with the electric field via the electric 
dipole moment p. The interaction Hamiltonian between 
the mesoscopic object and the optical mode oscillator 
then takes the form 



H lnt = -p-E oc (0) = -fQ. 



(43) 



so that the force on the oscillator coordinate is due to 
the electric dipole moment. The force is 

£Joo(n • p) = -y/ (47r/L 3 ) ujoop. (44) 

Eq.(44) follows from Eqs.(40), (42) and (43). 

If the oscillator were driven by an external force at 
frequency u>, then from Eqs.(25) and (26) the damping 
force and displacement at frequency u) obey 



8f u = n(w + iO + )5Q u 



(45) 



The electric dipole moment response to the electric field 
defines the polarizability a(C) of the mesoscopic object 
via 



S P u = a{oj + i0 + )n-<5E oow (0). 
Equivalently, 



Sp u = a(u + iO+)v/(47r/L 3 ) uj^SQ^ 



(46) 



(47) 



where Eqs.(40), (42) and (46) have been employed. From 
Eqs.(44), (45) and (47), one can relate the self energy 



part of the mode propagator to the polarizability of the 
mesoscopic object; The relationship is 

n(C) = «(0- (48) 

The condition for the oscillator stability 

loI = o£ - 11(0) > 0, (49) 

may be written in terms of the static polarizability 

47ra(0) < L 3 (stable oscillator). (50) 

For a very highly polarizable mesoscopic object, say a 
single domain ferroelectric grain [^l],[53 34 1, the photon 
oscillator can exhibit an instability. The general nature 
of the non-linear forces which eventually stabilize the os- 
cillator must come into play. 



V. SUPER-RADIANT PHASE 

The problem of how an oscillator reaches a new 
thermodynamic equilibrium position in the "unstable" 
regime can be discussed only in the context of a specific 
microscopic environmental Hamiltonian H which enters 
into Eq.(7). An often discussed model for a photon os- 
cillator employs N identical "two level atoms" p5|-|3l| in 
the thermal environment. The Hamiltonian is 



H 



N 



-eS z 



(51) 



where e is the energy difference between the single atom 
ground state and the excited state, <jj = {<Jxj,Oyj,&zj) 
are the Pauli matrices for the j th two level atom and 
S = (1/2) J2f=i a j 1S formally the "total spin" of the 
mesoscopic object. The object is placed in the region 
of the cavity where the electric fields of the photon os- 
cillator are concentrated. The electric field is presumed 
to interact with the atomic system via the total electric 
dipole moment component in the n direction; i.e. 



N 
3 = 1 



2fiS x , 



(52) 



where /i is the dipole moment matrix element. Using 
Eqs.(44) and (52) the force exerted on the photon oscil- 
lators by the dipole moment of the atomic system is 



16tt 

-jj- {uioofibx). 



(53) 



One is now in a position to compute the free energy 
G(Q,T) for all domains (stable or unstable) of the (now 
non-linear) oscillator. 
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Substituting Eqs (51) and (53) into Eq.(7) and taking 
the trace over the formal spin degrees of freedom yields 



G(Q,T)=uj 2 00 (Q 2 /2) + k B T In cosh 



2k B T 



k B T In cosh 




2k B T 



(54) 



Using Eqs (11) and (54), the shift in the frequency of 
the photon oscillators from Woo down to the temperature 
dependent ojq is found to be 



„ o(T) = ^ 1 _(?^V) tanJl (_i_). <„, 

The stability of the oscillator requires u)o(T) > 0, which 
in virtue of Eq.(55) reads 

T>T C (stable oscillator), (56) 

where the critical temperature for stability obeys 



tanh 



(■ 



L 3 e 



2k B T c ) \8irNu 2 



(57) 



The shift in the equilibrium position required to sta- 
bilize the oscillator when T < T c is found from the free 
energy minimum condition {dQ / dQ) T = 0. The resulting 
equilibrium position Q e obeys the transcendental equa- 
tion 



tanh 




1+ /lerro&^Qr 



L 3 e 2 



\ L 3 e 2 



(58) 



In Fig. 1, with y = (e/2k B T c ) = 0.1, we show 
the following: (i) Above the critical temperature x = 
(T/T c ) > 1 the renormalizcd frequency of the stable os- 
cillator uj (T) = v{x,y)to 00 is lowered as the temper- 
ature is lowered, (ii) Below the critical temperature 
x = (T/T c ) < 1, the oscillator position shift Q e (T) = 
(V Nk B T c /uj co )((x, y) (required to restore stability) in- 
creases as the temperature is lowered. Thus as the tem- 
perature is lowered from above the critical temperature, 
the oscillator frequency decreases until the "soft mode" 
has vanishing frequency at the critical temperature. The 
transition to an ordered phase begins as the temperature 
is lowered below the critical temperature. The oscillator 
seeks a new equilibrium position Q e {T < T c ) displaced 
from the old equilibrium position of Q = 0. The further 
the temperature is lowered, the more growth is exhibited 
in Q e . 

In Fig. 2 the free energy is plotted as a function of 
the displacement of the oscillator for two values of the 



temperature. For T > T c , the potential has a single 
minimum at Q = 0. For T < T c , the potential has a 
double minimum at Q — ±Q e (T). The oscillator will 
choose one of these minima in the super-radiant phase. 




x = (T/T c ) 

FIG. 1. For y = (e/fesT c ) = 0.2, the renormalized oscil- 
lator frequency luq = u) x v is shown above the critical tem- 
perature x = (T/Tc) > 1- The shifted equilibrium position 
Qe = (VNkBTc/uJoo)^ is plotted below the critical tempera- 
ture x = (T/T c ) < 1. 




FIG. 2. For y = (e/k B T c ) = 0.2, the free energy of the 
photon oscillator G(Q,T) = NksTg(^,x,y) is plotted as a 
function of Q = (V NksT /lo x )^ for two values of x = (T/T c ), 
above the critical temperature x > 1 and below the critical 
temperature x < 1. 



If a highly polarizable mesoscopic grain is coupled 
strongly to the transverse electromagnetic field via the 
electric dipole moment, then there is a transition from a 
normal phase (T > T c ) to a super radiant phase (T < T c ). 



VI. CONCLUSIONS 

It has been shown that any simple harmonic oscillator 
subject to a random force (producing damping) can be 
made unstable if the induced damping of the oscillator is 
sufficiently strong. If Woo is the frequency of the unper- 
turbed oscillator and if II(£) in Eq.(31) is the self energy 
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induced by the damping coefficient 3?er(w + i0 + ), then 
the precise condition for the instability is that 11(0) > ux^. 
It is possible to observe this instability for a mesoscopic 
grain which is highly polarizable. The grain must be 
placed in an optical cavity with a mode of spatial extent 
L large on the scale of the grain size. The condition for 
an optical cavity mode to go unstable is that 

47ra(0) > L? (unstable cavity mode), (59) 

where ol{oj + i0 + ) is the dynamic polarizability of the 
grain. The inequality in Eq.(59) can surely be reached 
with a single domain highly polarizable ferroelectric 
grain. 

The manner in which the oscillator is eventually sta- 
bilized (when Eq.(59) holds true) is not a trivial prob- 
lem. The example of a grain which consists of "two-level 
atoms" was worked out in detail. The normal and super- 
radiant free energies were exhibited in Fig. 2. When the 
frequency of the oscillator mode softens towards zero as 
T — > T c + + , the equilibrium position of the oscilla- 
tor shifts to a nonzero value Q e (T < T c ). For two level 
atoms, the details of the super-radiant transition can be 
computed. 

It is to be hoped that the theoretical considerations 
here put forth for ferroelectric grains will soon be put to 
laboratory tests. 
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